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ABSTRACT 

This paper describes the first optical spectroscopic survey of class I sources 
(also known as embedded sources and protostars) in the Taurus-Auriga dark 
cloud. We detect 10 of the 24 known class I sources in the cloud at 5500-9000 
A. All detected class I sources have strong Ha emission; most also have strong 
[O I] and [S II] emission. These data - together with high quality optical spectra 
of T Tauri stars in the Taurus- Auriga cloud - demonstrate that forbidden 
emission lines are stronger and more common in class I sources than in T Tauri 
stars. Our results also provide a clear discriminant in the frequency of forbidden 
line emission between weak-emission and classical T Tauri stars. In addition 
to strong emission lines, three class I sources have prominent TiO absorption 
bands. The M-type central stars of these sources mingle with optically visible 
T Tauri stars in the HR diagram and lie somewhat below both the birthline for 
spherical accretion and the deuterium burning sequence for disc accretion. 

Subject headings: stars: pre-main-sequence - stars: formation - ISM: jets and 
outflows 



1. INTRODUCTION 



Examining the earliest phases of low mass stellar evolution requires observations of 
protostars deeply embedded in the dense cores of nearby molecular clouds. These "class 
I" sources ( |Lada 1987| ) have blackbody-like spectral energy distributions that peak at 
wavelengths of 30-100 /im and bolometric luminosities, Lb ~ 0.1-100 Lq ( |Adams et al. 



[I987I ; llVlyers et al. 19S7[ [Wilking et al. 1989| ; [Kenyon et al. 1990| ; [Greene et al. 1994| ). 



Despite many detailed studies of their circumstellar environments (see, for example 
\et al. 1991a| ,b; [Andre et al. 1993| ; ['lerebey et al. 1993[ ; [Moriarty-Schieven et al. 1992 



amura 



1995; [Bontemps et al. 1996[ ; [Gomez et al. 1997[ ; [Hogerheijde et al. 1997| ; [Whitney et~al\ 



1997 ), understanding the stellar physics of these objects has proved elusive. Comparisons of 



observed bolometric luminosity functions with models is straightforward but controversial 
( [Wilking et al. 1989[ ; [Kenyon et al. 1990[ ; [1^'letcher fc Stahler 19944 1994b). The apparent 
lack of photospheric features in many objects has led several groups to abandon the HR 
diagram as a means for testing stellar evolutionary tracks of the youngest stars. These 
groups have proposed the bolometric temperature ([Myers et al. 19"98| and references 
therein), the submillimeter flux ([Saraceno et al. 1996[) , and the visual extinction ( [Adams 
[1990[) to replace effective temperature and have developed models to place evolving 
pre-main-sequence stars in their modified evolutionary diagrams. The accuracy of these 
techniques remains uncertain, because the methods are new and incompletely tested. 

In this paper, we report an optical spectroscopic survey designed to detect photospheric 
absorption features from the central stars of class I sources in the Taurus-Auriga cloud. 
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Although the hne-of-sight extinction to the central star is large, Ay ~ 30-60 mag ( [Whitney 



et at. 1997D , large ground-based telescopes can detect optical light scattered off cavities 
in the infalling envelopes of many objects. Optical data also provide the best measure of 
spectral types for pre-main-sequence stars. In general, I-band and J-band data are least 
contaminated by emission from an accretion disc and its associated boundary layer or 
magnetic accretion column (|Kenyon fc Hartmann 1990|). However, the very large continuum 



veiling detected on near-IR spectra of some class I sources ( [Uasali fc Matthews 1992 ; Greene 



fe Lada 1996a| , 1996b) favors I-band spectra, because the J-band veiling can be large if the 



disc extends to the stellar photosphere ([Kenyon fc Hartmann 19^ ; [Kenyon et al. 1996| ) 



Finally, optical spectra of class I sources allow an unambiguous comparison with optically 
brighter T Tauri stars, which have known spectral types in a well-calibrated system (see, 
for example, Kenyon fc Hartmann 1995; KH95 hereafter). 

Our results provide the first optical detection of M-type absorption features in an 
embedded protostar. We identify TiO absorption bands in three Taurus-Auriga class I 
sources; one other star may have TiO features and a fifth star may have K-type absorption 
features. We use optical spectra of T Tauri stars to calibrate the spectral types of class I 
sources and then construct a complete HR diagram for the Taurus-Auriga cloud. These 
data, coupled with new evolutionary tracks for protostars accreting from discs and two 
spectral types derived from near-IR spectra ( [Greene fc Lada 1996 b| ), show that class I 
sources in Taurus-Auriga mingle with T Tauri stars and lie below the birthline in the HR 
diagram. 

We also detect strong emission lines on the spectra of all protostars. Forbidden 
emission from [N II] and [S II] is much more common among class I sources than older, 
optically brighter stars having the same bolometric luminosity. The fluxes of forbidden 
emission lines also seem stronger among class I sources than other pre-main-sequence stars 
in the cloud. We find no evidence that the permitted emission lines, such as Ha and He I, 
are more common or stronger than in T Tauri stars. These results extend and confirm 
previous conclusions that jet activity declines as a pre-main-sequence star contracts to the 
main- sequence. 

We describe our observations in Sec. 2, explain our results in Sec. 3, and conclude with 
a brief discussion in Sec. 4. 



2. OBSERVATIONS 

We acquired optical spectra of faint Taurus- Auriga class I sources and other pre-main- 
sequence stars with the Red Channel Spectrograph at the Multiple Mirror Telescope (MMT; 
16-19 November 1995) and the Double Spectrograph at the Palomar 5-m telescope (29-30 
November 1995). At the MMT, we used the 270 g mm~^ grating and a 1" slit to produce 
spectra covering 5700-9000 A on a 1200 x 800 Loral CCD. On-chip binning of the pixels, 
2x2, yielded a spectral resolution of 10.8 A and a spatial resolution of 0''6 per pixel. At 
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Palomar, we used a 316 g mm~^ grating, a 2" slit, and a 1024 x 1024 CCD. The Palomar 
spectra cover 6000-8500 A with a spectral resolution of 10 A and a spatial resolution of 
0''47 per pixel. 

We obtained low resolution spectra of brighter young stars in Taurus-Auriga during 
1995-1996 with FAST, a high throughput, slit spectrograph mounted at the Fred L. 
Whipple Observatory 1.5-m telescope on Mount Hopkins, Arizona ([Fabricant et al. 1998| ). 



We used a 300 g mm ^ grating blazed at 4750 A, a 3" slit, and recorded the spectra on a 
thinned Loral 512 x 2688 CCD. These spectra cover 3800-7500 A at a resolution of ~ 6 A. 

We derive final object spectra using standard tasks within NOAO IRAF. After 
trimming the CCD frames at each end of the slit, we correct for the bias level, flat-field 
each frame, apply an illumination correction, and derive a full wavelength solution from 
calibration lamps acquired immediately after each exposure. The wavelength solution for 
each frame has a probable error of ±0.5-1.0 A. We extract object and sky spectra using 
the optimal extraction algorithm within APEXTRACT. We vary the size of the object/sky 
aperture from source to source to include additional radiation from extended emission. 

At both the MMT and Palomar, we selected class I sources based on published VRI 
photometry (see KH95), POSS red plates, and red narrow band continuum images ( Gomez 



\et al. 1997|) . Our sample of 10 class I sources is complete to R ~ 20: we detected all sources 
with optical counterparts on the POSS and several sources with measured R ~ 19-20. We 
did not observe 14 other class I sources without known optical counterparts and cannot 
estimate a reliable detection frequency for deeply embedded sources. The observed sample 
spans the observed range of class I luminosities in the Taurus-Auriga cloud. Lb ~ 0.2-20 
Lq, but does not include the reddest systems that do not have optical counterparts (see 
Whitney et al. 1997]). 



We selected FAST sources from the KH95 sample of ~ 150 known pre-main-sequence 
stars in the Taurus-Auriga cloud. We observed essentially all targets with V <^ 16-17: 
55 out of 65 weak emission T Tauri stars, 69 out of 96 classical T Tauri stars, and one 
class I source with a bright optical counterpart (Haro 6-28). Aside from their apparent 
brightness, these samples do not appear to be biased against any particular observational 
property of pre-main-sequence stars: the distributions of Ha equivalent widths, K-L colors, 
and bolometric luminosities of stars in the FAST sample are indistinguishable from the 
distributions for stars not included in the FAST sample using data from Cohen & Kuhi 
(1979) and KH95. We thus conclude that the FAST objects are a representative sample of 
known pre-main-sequence stars in the cloud. 



3. RESULTS 
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3.1. Basic Properties and Spectral Types 

Figure 1 shows FAST spectra of T Tauri stars with a range of emission characteristics. 
Weak emission T Tauri stars, such as LkCa 3, have absorption spectra similar to normal 
main-sequence stars with additional weak Ha emission lines ( |Uohen fc Kuhi 1979| ; |Walter 



at. 1987| ). Classical T Tauri stars, such as BP Tau, DG Tau, and DP Tau, have prominent 



emission lines and a variable blue continuum superimposed on a late-type absorption 
spectrum QCohen fc Kuhi 1979 ; Hartigan et al. 1989| , 1991). Some T Tauri stars have few 



emission lines other than H 1 and He I; others have prominent [O I], [S II], and [Fe II] 
emission lines ( Pohen fc Kuhi 1979| ). In most interpretations, the H I and He I lines form 
in an accretion region or outflowing wind close to the central star; the [O I] and [S II] lines 
form in a jet or in the wind (Pasri fc Bertout 1993| ; [Edwards et al. 1987| ; [Hamann fc Persson 
19921 ; see also [Martin 1997]) . 



Figure 2 shows contour maps of MMT spectra for several sources centered on Ha and 
[S II] AA6717,6730. We fit gaussian profiles to the continuum and a few emission lines 
along the spatial direction for each class I source and three point sources. Several class I 
sources - such as L1489 IRS (04016+2610), HH31 IRS2 (04248+2612), and 04264+2433 - 
are clearly extended along the slit, with a ~ 5"-10" compared to a ~ 1''0±0''2 for point 
sources. These sizes are comparable to sizes inferred from optical and near-infrared images 
of these sources ( [Whitney et al. 1997| ; [Gomez et al. 19Wi\ , and references therein). The 



very deeply embedded sources L1527 IRS (04368+2557) and L1551 IRS5 (04287+1801) are 
much more extended than other class I sources. Both have large optical reflection nebulae, 
^ 30" across, with multiple emission knots (see, for example, [Stocke et al. 198"^ ; [Graham~S 
Rubin 1992[ ; [Eiroa et al. 1991 [Gomez et al. 1997[) . 



Our results indicate that the continuum and emission lines of all class I sources are 
equally extended within the errors of the fit. The mean difference in spatial extent between 
the emission lines and the continuum is (cr) = 0''4 ± 0"3 for the 5 sources with strong 
continua (L1489 IRS, 04158+2805, HH31 IRS2, 04264+2433, and 04489+3042). This 
difference is small compared to the typical spatial extent of each source, a ~ 5"-10". 
Our observations have insufficient spatial resolution to discriminate between sources with 
emission knots such as L1527 IRS and those without emission knots such as 04489+3042. 

Figures 3-5 show MMT and Palomar spectra for Taurus class I sources. We detect 
a strong continuum in several objects (Figure 3), including L1489 IRS, HH31 IRS2, 
04264+2433, and 04489+3042. Three class I sources - 04158+2805, 04248+2612 and 
04489+3042 - have the deep TiO absorption bands characteristic of M-type stars. The 
continua of L1489 IRS and 04264+2433 appear featureless, although both have a prominent 
dip at 8100 A. All of these class I sources have a strong Ha emission line, along with 
moderately strong emission from [O I], [S II], and Ca II. Several have He I emission at 
AA5876, 6678, 7065. These spectra are similar to the spectra of classical T Tauri stars in 
Figure 1. 



The spectra of class I sources in Figure 5 more closely resemble spectra of jets or 
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Herbig-Haro (HH) objects (e.g., Bohm fc Soil ; Dopita et at. 1982 ; C^oodrich 198(: ; Reipurtli 



fc Heathcote 1991| ; see also |Raga et al. 1996| ). We detect little, if any, continuum emission 



from 04239+2436, 04295+2251, and L1527 IRS, but all of these objects have very strong 
emission lines (Figure 5). The relative intensities of the emission lines in 04239+2436 
and L1527 IRS are similar to those for on-source emission in HH30 IRS (Figure 4; right 
panel) and L1551 IRS5 (Figure 5; upper right panel). We detect [O I] and [S II] emission 
in all class I sources except 04295+2251 (Figure 5; lower left panel), where we identify a 
prominent Ha emission line superimposed on a very weak continuum. 

To compare spectra of class I sources and T Tauri stars in more detail, we measure 
the strengths of several prominent absorption and emission lines. We fit gaussian profiles 
to obvious emission lines using SPLOT within NOAO IRAF and use the deblend option for 
blended lines such as the [S II] doublet and the He I A5876 and Na I blend. In the absence 
of accurate de-reddened fluxes, the equivalent width, EW, provides a good relative measure 
of emission line strengths for strong continuum sources. We place upper limits of 100-200 
A for equivalent widths of weak continuum sources depending on the continuum level. 
Table 1 lists our results for class I sources; Table 2 summarizes measurements for FAST 
spectra of T Tauri stars. We include separate entries for the MMT and Palomar spectra of 
two objects, L1489 IRS and HH31 IRS2, and list the MMT results first in both cases. We 
estimate probable errors of ±10% for strong lines with EW > 10 A and ±20% for weaker 
lines based on measurements of 2-3 separate spectra for each star in the FAST sample and 
two stars in the MMT/Palomar sample. 

We measure spectral types using TiO absorption indices, defined as the depth of a 
TiO band at a wavelength A relative to an interpolated continuum point ( P'Connell 1973| ). 
The TiO bands at A6180 and A7100 are temperature-sensitive for M dwarfs and giants 
( P'Connell 1973|) . Kenyon & Fernandez- Castro (1987) derive reliable spectral types for 
the red giant components in symbiotic stars - which also have strong emission lines and a 
variable blue continuum - with these features. We define: 



[TiO], = -2.5 log [ ,n22^7r F^) 

V-T 6125 + U.Z/:£)(^i:' 6370 ~ -^6125^/ 

and 

[TiO], = -2.5 log ( J;^°° —-\ . (2) 

V-T 7025 + U./l^r 7400 — -r 7025j/ 

Each 30 A bandpass used for these indices avoids contamination from strong emission lines 
and telluric absorption bands. The two TiO indices increase from [TiO]i ~ [TiO] 2 ~ at 
K4-K5 spectral types to [TiO]i ~ 0.8 and [TiO] 2 ~ 1.0 at M6 spectral types. 

Figure 6 shows [TiO] 2 as a function of [TiO]i for normal main-sequence stars (filled 
circles), T Tauri stars with weak emission lines on FAST spectra (light triangles), and 



- 7- 



pre-main sequence stars with MMT or Palomar spectra (crosses and diamonds). The locus 
of T Tauri stars generally follows the main-sequence stars except near [TiO]i ~ 0.5, where 
pre-main- sequence stars have larger [TiO]2 indices compared to main-sequence stars. Most 
pre-main- sequence stars with MMT or Palomar spectra lie on the T Tauri star locus. 
Both L1489 IRS and 04303+2240 have featureless continua and negligible TiO absorption 
on their spectra. Other stars with strong optical continua have modest to strong TiO 
absorption bands and must have M-type central stars. 

The measured TiO indices indicate optical veiling in two sources with MMT or Palomar 
spectra. The class I sources HH32 IRS2 and 04158+2805 lie above the pre-main-sequence 
locus in Figure 6. In both cases, the M-type absorption features for A > 7400 A are 
very strong, which suggests that the [TiO]i index is 'weak' compared to the [TiO]2 index. 
Many T Tauri stars with strong emission lines also have weak [TiO]i indices. These T 
Tauri stars have substantial emission from a blue continuum source which veils optical 
absorption lines; this veiling increases towards short wavelengths in all cases (see [Hartigan 
\et al. 19911) . Optical veiling from a hot, T ~ 10'^ K, continuum source probably causes the 
weak [TiO]i index in HH32 IRS2 and 04158+2805, but our optical spectra have insufficient 
signal-to-noise to verify that absorption features at A < 6000 A are similarly weakened. 

We derive spectral types for class I sources and faint T Tauri stars via comparison with 
T Tauri stars of known spectral type. We adopt spectral types for bright T Tauri stars 
from KH95 and use FAST spectra to calibrate [TiO]i and [TiO] 2 as a function of spectral 
type. The measured TiO indices for class I sources then yield the spectral types listed in 
Table 1. We estimate probable errors of ±1-2 subclasses for the spectral types based on 
measurement errors of the TiO indices and the intrinsic uncertainty in assigning spectral 
types to bright T Tauri stars. 



3.2. Jet Emission 

As we noted in the introduction, powerful optical and molecular outflows distinguish 
class I sources from other pre-main-sequence stars in nearby dark clouds. Practically 
all class I sources have molecular outflows; very few optically visible T Tauri stars are 
associated with molecular outflows (e.g., [Bontemps et al. 1996|) . In their complete survey 
of Taurus class I sources, Gomez et al. (1997) show that the frequency of optical jets 
decreases from ^ 60% for class I sources to <^ 10% for T Tauri stars. Among T Tauri stars, 
optical jet emission is almost always associated with classical T Tauri stars (CTTS) instead 
of weak-emission T Tauri stars (WTTS; see lildwards et al. 1987 ; Hartigan et al. 1995|) . 



Near-infrared surveys also indicate more emission line activity among class I sources than 
CTTS or WTTS ( [Greene fc Lada 19964 1996b). These results suggest that jet emission is 
correlated with disc accretion and that disc accretion somehow declines from class I sources 
to T Tauri stars ( [Edwards et al. 1987| ; [Hartigan et al. 1995[ ; preene fc Lada 1996b| ). 



To examine these correlations with our spectroscopic data, we divide sources into 
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classes based on the ratio of far- infrared to bolometric luminosity, Lpjji/L\^ (KH95). In this 
system, class I sources have Lpjji/L]^ > 0.8, CTTS have Lpir/L\^ ^ 0.1-0.3, and WTTS 
sources have Lpj^/L^, < 0.1. We also select 14 sources in Tables 1-2 with Lpui/Li, ^ 
0.3-0.8 as fiat-spectrum sources (see also KH95; Preene fc Lada 1996b| , 1997). All of these 



sources have the same median luminosity. Lb ~ 0.5-0.8 L0, except for the fiat-spectrum 
which have Lb ~ 1.5 L0 (KH95). For each class, we compute the detection frequency 
for each of the emission lines listed in Tables 1-2. Figure 7 shows our results for [S 11]. 
The frequency of [S II] emission obviously increases with increasing Lfir/L-^,- We find a 
similarly strilcing trend for [N II] emission: practically all class I sources have strong [N II] 
but only a few CTTS or WTTS have any [N II] emission. 

Our detection frequencies for forbidden-line emission are lower limits, because we 
cannot detect weak emission lines with EW <^ 0.2-0.4 A (see Tables 1-2). To estimate 
the importance of this uncertainty, we compare our results with Hartigan et al. (1995) 
who derive emission line equivalent widths from echelle spectra. At higher resolution, the 
detection frequency for [S II] emission increases to ~ 50% for CTTS and remains unchanged 
at 0% for WTTS. The [N II] emission shows a similar trend but is detected less often than 
[S II]. Hartigan et al. note, however, that the forbidden emission lines in CTTS consist 
of high velocity material from the jet and low velocity gas near the disc. The wealc, low 
velocity emission is responsible for the larger [S II] detection frequency among CTTS in the 
Hartigan et al. sample. Our low resolution spectra do not detect this emission and thus 
provide a proper estimate for the frequency of high velocity jet emission among CTTS. 

To check further the reality of the trend in Figure 7, we perform a simple test. We 
assume a parent population of N sources with an intrinsic probability, p, of [S II] emission. 
For Uj observed sources, the probability of detecting [S II] emission in kj sources is given 
by the binomial distribution: 



kj\ [rij -kj)] 

for rij < N. If we require Pobs ^ 10^^, the allowed range in p for a single class of 
pre-main-sequence star is large: p <^ 0.13 for WTTS, p ~ 0.03-0.36 for CTTS, p ~ 
0.24-0.93 for fiat-spectrum sources, and p ~ 0.33-1.00 for class I sources. However, the 
probability of realizing the observed frequency of [S II] detection for any two classes from a 
single parent population is extremely small. The probability that the WTTS and CTTS in 
our sample have the same parent distribution never exceeds 10"'^; it exceeds 10~^ only for 
p = 0.04-0.12. We find no common intrinsic probability for the CTTS and fiat spectrum 
sources or for CTTS and class I sources: Pots is less than 10~^ for any value of p. The 
observed frequencies of [S II] emission in class I and fiat spectrum sources, however, could 
be chosen from the same parent population for p = 0.46-0.90 if pobs > 10~^. 

We derive similar results for the frequency of [N II] emission. We detect the AA6548, 
6584 doublet in 82% of 11 class I sources, 36% of 14 fiat spectrum sources, 4% of 46 CTTS, 
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and 0% of 54 WTTS with reliable Lpir/L\^. The allowed ranges in the intrinsic probability 
of [N II] emission are then p ^ 0.33 for class I sources, p ~ 0.07-0.76 for flat-spectrum 
sources, p ^ 0.21 for CTTS, and p ^ 0.12 for WTTS. The observed detection frequencies 
allow a single intrinsic probability for [N II] emission between CTTS and WTTS for p ^ 
0.11 and between class I and flat spectrum sources for p ^ 0.40-0.71. Our data do not allow 
a single [N II] emission probability for all sources. This result is not as strong as for the 
[S II] lines, because weak [N II] emission is more difficult to detect due to the strong Ha 
lines in many CTTS. 

We conclude that the increasing frequency of forbidden line emission as a function 
of Lfir/L\^ is real. Sources with large Lpf^/L]^ are much more likely to be associated 
with forbidden line emission than sources with small Lpjji/L]^. In addition, our sample is 
large enough to detect a significant difference in the frequency of [S II] emission between 
CTTS and WTTS. Previously published data had suggested this difference, but the data 
were too heterogeneous to make a firm statistical comparison. Our sample is not large 
enough to detect a difference in the [N II] or [S II] emission frequency between class I and 
flat spectrum sources. Adding sources to the optical sample and improved measurements 
of Lpifi/L]^ would allow a better discriminant between class I and flat spectrum sources. 
These improvements require larger ground-based telescopes and new far-IR data from either 
ISO or SIRTF. 



3.3. The HR Diagram 

Figure 8 shows an HR diagram for Taurus-Auriga pre-main- sequence stars (see 
KH95). Filled circles plot WTTS and CTTS from KH95Q. Crosses indicate the positions 
of 04158+2805, HH31 IRS2, and 04489+3042 using our new spectral types. The triangle 
indicates Haro 6-28 with a revised spectral type (M2) based on FAST spectra. The 
diamonds denote 04181+2655 and 04295+2251 using near-IR spectra from Greene & Lada 
(1996b). The relative positions of these two class I sources should be accepted with some 
caution, because near-IR spectral types for 04489+3042 (K1-K2) and Haro 6-28 (K5-K6) 
are much earlier than our optical spectral types. 

Although the class I sample is small, the positions of class I sources are not especially 
distinct from the distribution of WTTS and CTTS. Three class I sources straddle the 10^ 
yr isochrone from the CMA models of D'Antona & Mazzitelli (1994). Two other class I 
sources fall midway between the 10^ yr and 10^ yr isochrones. All class I sources are within 
2a of the 10^ yr isochrone and the stellar "birthline" for spherical accretion (Stahler 1983, 



1988; IPalla fc Stahler 199q , 1993; see also KH95). Only Haro 6-28, however, lies on the 



We use the KH95 conversion from spectral type to effective temperature for the class I sources. 
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stellar birthline^ 

To compare our new data with an alternative to spherical protostellar accretion theory 
( [Stahler 198ij , 1988), we consider the evolution of protostars accreting material from a disc. 



We construct a set of stellar models using the most recent version of the Eggleton evolution 
program ([Eggleton 1971| , 1972, 1973). Our models assume an initially uniform composition 



with abundances of hydrogen X = 0.7, helium Y = 0.28, deuterium Xd = 3.5 x 10~^, and 
metals Z = 0.02 appropriate for the meteoritic mixture determined by Anders & Grevesse 
(1989). Pols et al. (1995) describe the equation of state, which includes molecular hydrogen, 
pressure ionization, and coulomb interactions. The nuclear reaction network includes the 
pp chain and the CNO cycles. Deuterium burning is explicitly included at temperatures too 
low for the pp chain. Once the pp chain is active hydrogen burns to He^ with deuterium 
and He^ in equilibrium. The burning of He'^ is not explicitly followed. We use the opacity 
tables of Iglesias, Rogers & Wilson (1992) and Alexander & Ferguson (1994). We adopt an 
Eddington approximation ( [Woolley fc Stibbs 1953|) for the surface boundary conditions at 
an optical depth of r = 2/3. Low-temperature atmospheres, in which convection extends 
out as far as r ^ 0.01 ( [Baratt'e et at. 19"95| ), are not modeled completely. However the effect 



on observable quantities is not significant (see [Kroupa fc 'lout 1998| ). 



In these calculations, the initial protostar is a fully convective 0.1 object with a 
radius of 3 Rq and an effective temperature of lO'^'^^ K. This starting point lies just off 
the right boundary of Figure 8. We add accreted material to the surface with the initial 
composition and with the same state (entropy, temperature, etc.) as the surface. In this 
approximation, most of the stellar surface is free to radiate normally with the boundary 
conditions described above. These boundary conditions are a compromise, because we are 
modeling a two dimensional process with a one dimensional evolution code. 

The solid lines in Figure 8 indicate two accreting protostellar tracks using the Eggleton 
code. The thinner line accretes at lO~^M0yr~^; the thick line accretes at 10~^MQyr~^. 
Initially, both stars contract at fairly constant luminosity and move to the left in the HR 
diagram. This contraction continues until deuterium ignites at their centers. The tracks 
then turn upwards to follow a deuterium burning sequence. Both protostars are fully 
convective, so newly accreted material replenishes central deuterium. These protostars 
thus remain on the deuterium burning track until either accretion ceases or the rate of 
deuterium replenishment becomes insufficient for burning to continue to support the star. 
The latter occurs first, temporally, for the 10~^ yr~^ model at Tcfj = 10^'^^ K when it has 
reached a mass of 0.686 M©. The lO~^M0yr~^ model drops below the deuterium sequence 
at Tgff = lO'^-^^ K at a mass of 0.456 Mq. In both cases, we continue accretion until the total 
mass reaches 0.96 Mq. During the protostellar evolution, the stars lie close to the Hayashi 
tracks appropriate to their instantaneous mass. If accretion ceases at any time, the star will 



^We suspect that Haro 6-28 may be a close binary similar to GV Tau (Haro 6-10), because their spectral 
energy distributions arc similar. Both GV Tau and Haro 6-28 are 2-3 mag brighter in the optical than other 
class I sources in the sample. 
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shrink down to the main- sequence along a normal pre-main-sequence track and will ignite 
and burn deuterium on the way if it has not already done so. 

In our accretion models, the deuterium burning sequence defines a stellar birthline in 
the HR diagram. Previous calculations find similar results. Stabler (1983, 1988; see also 
Palla fc Stabler 1990| , 1993) first identified the birthline for spherical accretion and showed 



that this locus provides a good upper envelope for observations of young stars in nearby 
molecular clouds. Mercer-Smith et al. (1984) published the first HR diagram track for 
a disc-accreting protostar using a code similar in spirit to the Eggleton code. Hartmann 
et al. (1997) later derived a birthline for disc accretion from semi-analytic calculations. 
Our deuterium burning sequence lies close to Stahler's (1983, 1988) birthline for log Tgg ^ 
3.55 and falls 5 log L k, 0.1-0.2 below the birthline for log Teg ^ 3.55. The displacement 
refiects differences in the starting conditions and outer boundary condition. The birthlines 
converge at large Tgff, because the boundary conditions become less important as the stellar 
luminosity increases (see also Palla fc Stabler 1990| , 1993). 



The location of the stellar birthline in our models is sensitive to the adopted deuterium 
abundance. This behavior follows from the explicit dependence of the stellar luminosity 
on the rate of deuterium burning (see, for example, [Stabler 1988| ; [Hartmann et al. 1997| ). 



Our deuterium burning sequence shifts by 5 log L —0.13 for a factor of two reduction 
in Xj). Our birthline then roughly coincides with that of Hartmann et al. (1997) for 
Xd = 1.75 X 10-5. 

With only six class I sources in our HR diagram, it is difficult for the data to favor 
convincingly any theoretical calculation. The 2o" error bars are consistent with all of the 
tracks, even without considering uncertainties in the model input parameters. The data 
lie closer to the disc-accretion birthline of Hartmann et al. (1997) and our deuterium 
burning sequence than either the birthline for spherical accretion or the accretion track of 
Mercer-Smith et al. (1984). In all cases, changing model input parameters - such as the 
deuterium abundance - would allow a better match between data and the models. More 
rigorous comparisons thus await observational estimates of unknown quantities such as the 
deuterium abundance and a larger sample of class I sources with rehable spectral types and 
luminosities. 

We conclude this section with several points about the comparison of observations 
with model tracks. First, the distribution of class I sources about any birthline should be 
uniform, if the range in initial conditions is small. This dispersion should be comparable to 
the observational errors. In our case, 5 out of 6 class I sources fall l-2cr below the birthline 
for spherical accretion ( [Stabler 1983[ 1988) and the Mercer-Smith et al. (1984) accretion 



track. The data are somewhat more consistent with our deuterium burning sequence and 
the Hartmann et al. (1997) birthline for an accretion rate of 10"^ Moyr"^. A larger sample 
of class I sources should distinguish between models. 

Second, it is important to compare the stellar component of the protostellar luminosity 
with the predictions of model tracks. The observed class I luminosities are the total 
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luminosity and have not been corrected for the unknown amount of accretion luminosity. 
In most of the CTTS shown in Figure 8, the accretion luminosity is a small fraction - 
10% to 30% - of the stellar luminosity plotted in the figure (see [Hartigan et al. 1991| , 
1995; [GuUbring et at. 1998|). However, the accretion luminosity is roughly comparable to 



the stellar luminosity in the continuum + emission sources and is ~ 100 times the stellar 
luminosity in FU Ori systems like L1551 1RS5 ( [Hartigan et al. 1991| ; [Hartmann fc Kenyon 



1996 ). We expect a small accretion contribution for class I sources with optical spectra 



similar to most CTTS, although the accretion luminosity in L1489 IRS may be large (see 



also [Greene fc Lada 1996b| ). The arrow in Fig. 8 indicates the displacement in log T^s for 
a 5 log L = —0.1 change in the bolometric luminosity. This change moves the data closer 
to both the birthline for spherical accretion and our deuterium sequence and illustrates the 
difficulty in comparing models with current data. 

Third, knowledge of the deuterium abundance, and to a lesser extent the lithium 
abundance, is also necessary to compare observations with model predictions. The mean 
deuterium abundance for stars in a molecular cloud sets the location of the birthline in the 
HR diagram; any star-to-star scatter in the abundance spreads the birthline vertically in the 
HR diagram. To our knowledge, the deuterium abundance has not been measured in any 
pre-main-sequence star. Recent measurements indicate a factor of 3-5 scatter in the lithium 
abundance among nearby molecular clouds (e.g., [King 1995| ; [Lee et al. 1994[ ; puncan ^ 



RebuU 1996[ ) that could be due actual abundance differences between stars (see [Lee et al 



[1994[ ) or differences in the analysis procedures (see puncan fc RebuU 1996[ and references 



therein). Observations of older open clusters may also indicate considerable star-to-star 
differences in the rate of lithium depletion among stars with the same mass (e.g., [Randich 



\et al. 19971 and references therein). Similar spreads in the deuterium burning rate, due 



perhaps to star-to-star variations in accretion rate, further complicates the comparison of 
observations with model tracks. The sample of protostars is not currently large enough to 
worry about abundance variations, but the uncertainties will become more important as 
sample sizes increase. 



4. DISCUSSION and SUMMARY 



In the previous sections, we have described the first optical spectroscopic survey of 
class I, embedded sources in a single molecular cloud. We supplemented these data with 
high quality optical spectra of a representative sample of older and optically brighter T 
Tauri stars. The combined set of spectra show that the optical spectra of class I sources 
qualitatively resemble the optical spectra of T Tauri stars. Our analysis further reveals 
common physical properties and substantial differences between class I sources and T Tauri 
stars, as summarized below. 

Our data provide the first indication that the distribution of stellar spectral types 
among class I sources may not be very different from that of WTTS and CTTS. Of the five 
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class I sources with strong optical continua, one (L1489 IRS) is a continuum + emission 
source, three are M-type stars, and another (04264+2433) may have an M-type central 
star. To the best of our knowledge, these are the first low mass protostars with measured 
optical spectral types. This sample is too small for a meaningful comparison with the 
distribution of spectral types among more evolved pre-main-sequence stars in the cloud. 
We note, however, that the median spectral type for WTTS and CTTS is K7-M0 and that 
the frequency of continuum + emission sources is ~ 5%-10% (KH95). 

Published observations indicate other similarities between class I sources and older 
pre-main-sequence stars in Taurus-Auriga. First, class I sources have the same intrinsic 
near-IR colors as do CTTS. Whitney et al. (1997) show that the observed near-IR colors 
of class I sources can be modeled as a CTTS surrounded by an infalling envelope with 
an optical extinction. Ay ~ 30-60 mag. This analysis leads to the conclusion that the 
radiation from the star and inner disc of a class I source is similar to that of a T Tauri 
star (see also [Greene fc Lada 1996b| ; Calvet et al. 1997 reach a different conclusion). 



Second, the bolometric luminosity distributions of class I sources, CTTS, and WTTS 
are indistinguishable (KH95). All three groups of pre-main-sequence stars have median 
luminosities of Lb ~ 0.5-0.8 Lq. This unexpected result is supported by the positions of 
class I sources in the HR diagram. Our data show that class I sources have luminosities 
and effective temperatures very similar to those of CTTS and WTTS in the cloud. These 
conclusions are surprising, because a class I source should have a larger luminosity once it 
has accreted nearly all of its final mass, and this luminosity should decline with time as the 
star approaches the main sequence (see, for example, [Stahler 1983| , 1988; [Falla fc Stahler 



1993| ; [Hartmann et al. 1997| ; Fig. 8). The current sample, however, is too small to test 



stellar models in detail. The errors in luminosity and effective temperature are also too 
large. Observations with the next generation of ground-based telescopes will undoubtedly 
expand the sample, reduce the errors, and provide better tests of protostellar accretion 
theory. 

One feature that distinguishes class I sources is their strong forbidden-line emission. 
As a group, class I sources are much more likely to have forbidden-line emission than CTTS 
or WTTS (Fig. 7). This result confirms previous conclusions from imaging data (e.g., 
Gomez et al. 19971) and indicates that class I sources are more likely to drive outflows 



than CTTS or WTTS (see also [Bontemps et al. 1996|; [Moriarty-Schieven et al. 1992], 



1994). The equivalent widths of the forbidden lines are also larger in class I sources than in 
CTTS or WTTS. Although some large equivalent widths may be due to very weak optical 
continua, the [S II] equivalent widths in HH31 IRS2 - a class I source with a prominent TiO 
absorption band - are larger than observed in any CTTS in our sample (see Tables 1-2). 
Deeper optical spectra of our sample and other class I sources would clarify this point. 

Our sample is not large enough to test whether class I sources also have more prominent 
permitted emission lines than CTTS. The median Ha equivalent width for class I sources, 
~ 90 A, is much larger than the median equivalent width for CTTS, ~ 30-40 A. This 
difference is roughly what we expect if class I sources have larger continuum veiling than 
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CTTS ( lUasali fc Eiroa 1996| ; Greene fc Lada 1997|) and if the Ha equivalent width correlates 



with veiling ( [Hartigan et al. 1995| and references therein). However, the frequency of He I 



AA5876, 6678 emission among class 1 sources is roughly comparable to that among CTTS. 
We measure a He I emission frequency of 50% among 6 class I sources with reasonable 
signal-to-noise at 6000 A, 57% among 14 flat-spectrum sources, and 65% among 46 class 
II sources. For both emission lines, the class I sample is probably biased against small 
equivalent widths, because class I sources without emission lines are probably fainter than 
sources with emission lines. A deeper survey with a larger telescope could enlarge the 
sample of class I sources with high quality optical spectra. These data would provide a 
good test for differences in the distribution of Ha equivalent widths between class I sources 
and CTTS. 

These results fit into the general picture of Taurus-Auriga class I sources developed 
in KH95 and in Kenyon et al. (1990). In this picture, class I sources are envelopes of gas 
and dust falling into the central star-disc system at rates of a few x 10~^ Mq yr~^ (see also 
Adams et al. 1987| ; |Kenyon et al. 19934 1993b; [Whitney et al. 1997|) . Bell & Lin (1994) 



show that the stable accretion rate through the disc onto the central star is either very 
low - ^ a few X 10^'^ Mq yr^-*^ - or very high - a few x 10^^ yr"-*^ - compared to 
the infall rate. The disc spends most of its time in the low accretion rate state; the disc 
mass then slowly increases with time until it reaches a critical level and evolves to the high 
accretion rate state. This model explains the low observed luminosities of nearly all class I 
sources as well as the occasional high luminosity of a source such as L1551 IRS5. 

Models with time-dependent disc accretion also qualitatively account for the evolution 
of forbidden and permitted emission lines in pre-main- sequence stars. We expect the 
time-averaged accretion rate through the disc to decline as the envelope disperses. If the 
Ha and other permitted emission lines of class I sources form in the accretion region of the 
inner disc as in CTTS, then the median Ha equivalent width should decline as a pre-main 
sequence star evolves from a class I source into a CTTS and then into a WTTS. Most 
models for jet formation link the mass loss rate in the jet to the mass accretion rate through 
the disc (see, for example, Pabrit et al. 1990| ; [Najita fc Shu 199^ ; |Shu et al. 19944 1994b), 



so we expect forbidden emission to decline as well. Explaining the observations of emission 
line equivalent widths with a quantitative model of a dispersing envelope and evolving disc, 
however, is currently beyond our reach. 

Finally, our results further demonstrate the advantages of optical spectra. Recent 
surveys of larger samples of class I sources using near-IR spectroscopy have yielded only two 
spectral types each in Taurus-Auriga (|Casali fc Eiroa 1996| ; [Greene fc Lada 1996b| ) and p 



Oph ([Greene fc Lada 1996b| , 1997). Casah fc Eiroa (1996; see also [Casali fc Matthews 1992 
[Greene fc Lada 1996b| , 1997) conclude that continuum emission from dust in a circumstellar 
disc or envelope veils photospheric absorption features on near-IR spectra of class I sources. 
Preliminary results further suggest that this veiling is larger in class I sources than in CTTS 
or WTTS ( Pasali fc Eiroa 1996| ; [Greene fc Lada 1996b[ , 1997). Dust emission is much 



weaker relative to a normal stellar photosphere at shorter wavelengths, ^ 1 /im, so optical 
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spectra may yet provide the best measure of spectral types in class I sources. 
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Table 1. Spectral types and equivalent widths from MMT and Palomar spectra 



Source Name 


SED 


r / T 
^FIR/^h 


ST 


[O IJ 6300 


[O IJ 6363 


Ha 


Fat ttT f o a 

[N IIJ 6584 


rci TTl ^7 -1 n 

[S IIJ 6717 


[S IIJ 6730 


L1489 IRS 


I 


0.85 


Ce 






55.5 


1.5 


2.0 




L1489 IRS 


I 


0.85 


Ce 






56.0 


1.4 


1.3 


1.4 


04108+2805 


II 




MO 






37.0 








04113+2758 


II 


<0.77 


M2 






49.0 








CoKu Tau/1 


II 




MO 


71.5 


23.0 


110.5 


33.5 


20.0 


31.0 


04158+2805 


I 


<0.93 


M3 


91.0 


35.5 


326.0 


24.4 


13.5 


18.0 


04216+2603 


II 


<0.48 


MO 






26.5 








04239+2436 


I 


0.91 








>200.0 


>100.0 


>100 


>200.0 


HH31 IRS2 


I 


0.92 


M2 






86.0 


13.0 


8.0 


10.0 




T 
1 


n no 
0.9/ 


A/ro 
MZ 






95. U 


l.b 


D.O 


9.0 


04260+2642 


I 


0.89 




121.0 


38.0 


125.0 


14.5 


13.5 


24.5 


04264+2433 


I 


0.95 


MO? 


23.0 


7.0 


78.5 


10.5 


7.0 


10.0 


L1551 IRS5 


I 


0.98 


G-K 


41.0 


13.0 


83.0 


32.0 


40.0 


52.0 


LkHa 358 


II 


<0.76 


M5-M6 






86.5 








HH30 IRS 


I/II 




MO? 


110.0 


35.0 


185.0 


30.5 


47.0 


77.5 


04295+2251 


I 


0.80 








65.5 








04303+2240 


I 


<0.18 


Ce 


20.0 


7.5 


122.0 


16.5 






04368+2557 


I 


1.00 








>200.0 


>100.0 


>100 


>200.0 


04385+2550 


II 


0.72 


MO 






20.0 


3.0 


2.5 


3.5 


04489+3042 


I 


0.83 


M2 


164.0 


22.5 


11.0 


24.5 


7.0 


11.0 



Notes for Table 1: 

1. Three class I sources have detectable He I AA5876, 6678 emission with EW « 2-3 A: L1489 IRS, HH31 IRS2, and 
04264+2433. Two class II sources have comparably strong He I emission: 04216+2603 and 04303+2240. 

2. L1489 IRS is also known as 04016+2610. 

3. 04185+2805 is class II in KH95. Our optical spectra - together with unpublished optical and near-IR photometry - 
indicate that this object is a class I source. 

4. L1551 IRS5 is also known as 04287+1801. 

5. HH31 IRS2 is also known as 04248+2612. 
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Xaule Z 


Emission 


line equivalent widths from FAST spectra 


Source Name 


SED 


LfirIL\, 


Ha [N II] 6584 He I 6678 [S II] 6717+6730 


HBC 351 


HI 


0.05 


2.8 


HBC 352 


HI 


<0.03 


-2.2 


HBC 353 


HI 


<0.04 


-2.1 


HBC 354 


HI 


0.14 


-1.2 


HBC 355 


HI 


0.07 


-1.2 


HBC 356 


HI 


0.06 


0.8 


HBC 357 


HI 


0.10 


1.1 


HBC 358 


III 


0.04 


10.0 


HBC 359 


III 


0.04 


5.4 


HBC 360 


III 


<0.20 


6.7 


HBC 361 


HI 


<0.20 


6.4 


HBC 362 


HI 


<0.15 


6.6 


LkCa 1 


HI 


<0.05 


3.4 


Anon 1 


HI 


<0.02 


2.5 


04108+2910 


II 


<0.43 


37.5 ••• ••• 2.0 


V773 Tau 


II 


0.10 


2.8 


FM Tau 


II 


0.19 


101.0 ■■■ 1.6 


FN Tau 


II 


0.31 


16.0 


CW Tau 


II 


0.23 


140.0 ■■■ 1.2 1.1 


CIDA 1 


II 


<0.18 


149.0 


LkCa 3 


HI 


0.01 


2.7 


FP Tau 


HI 


0.13 


27.0 0.7 


FO Tau 


II 


0.15 


116.0 ■■■ 1.5 


CX Tau 


II 


0.17 


18.0 


CIDA 2 


HI 


<0.09 


6.8 


LkCa 4 


HI 


<0.02 


3.2 


CY Tau 


II 


0.10 


63.0 ■■■ 1.7 


LkCa 5 


HI 


0.16 


3.8 


CIDA 3 


II 


<0.14 


10.5 


V410 Tau 


HI 


0.01 


2.3 


DD Tau 


II 


0.40 


177.5 ••• 3.9 2.2 


CZ Tau 


II 


0.29 


6.7 


HBC 372 


HI 


<0.15 


-1.2 


V892 Tau 


II 


0.58 


6.6 


Hubble 4 
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FIGURE CAPTIONS 



Fig. 1. — Optical spectra of T Tauri stars. The degree of excitation increases clockwise from 
the weak emission line T Tauri star, LkCa 3, through the classical T Tauri stars, BP Tau 
and DP Tau, to the continuum + emission star DG Tau. 

Fig. 2. — Spatially resolved spectra of Taurus pre-main-sequence stars. The left panels plot 
Ha spectra for two class I sources, L1489 IRS (04016+2610) and HH31 IRS2 (04248+2612). 
The right panels plot [S II] spectra for L1527 IRS and HH30 IRS. All sources are very 
extended compared to the l"-2" extent of an unresolved point source. 

Fig. 3. — Optical spectra of four class I sources in Taurus-Auriga. All class I sources have 
strong emission from Ha, [O I], and [S II]; some also have He I and Ca II triplet emission. Two 
systems - HH31 IRS2 and 04489+3042 - have the deep TiO absorption bands characteristic 
of M-type stars. A telluric absorption feature at A7650 is present in the spectra of L1489 IRS 
and 04264+2433; this feature is blended with a TiO band in HH31 IRS2 and 04489+3042. 

Fig. 4. — Optical spectra of two Taurus-Auriga class I sources. One class I source, 
04158+2805, has deep TiO absorption bands and strong Ha, [O I], and [S II] emission 
lines. The spectrum of HH30 IRS is similar to other Herbig-Haro objects, with very intense 
emission from Ha, [O I], [N II], and [S II]. A telluric absorption feature at A7650 is present 
in the spectrum of HH30 IRS; this feature is blended with a TiO band in 04158+2805. 

Fig. 5. — Optical spectra of four class I sources in Taurus-Auriga. These objects have 
negligible continuum emission and very strong emission from Ha, [O I], [N II], and [S II]. 
The emission line equivalent widths are comparable to those observed in Herbig-Haro objects. 

Fig. 6. — TiO indices for field dwarfs and Taurus- Auriga pre-main-sequence stars. Several 
class I sources (diamonds) and heavily reddened T Tauri stars (plusses) have TiO band 
strengths comparable to those observed in T Tauri stars with negligible optical veiling (light 
filled triangles) and normal main-sequence stars (filled circles). 

Fig. 7. — Frequency of [S II] emission among Taurus-Auriga pre-main-sequence stars. The 
[S II] emission frequency increases with the ratio of far-IR to bolometric luminosity, Lpm/L^j. 

Fig. 8. — HR diagram for Taurus-Auriga pre-main-sequence stars. The filled circles are 
T Tauri stars from Kenyon & Hartmann (1995). The crosses are class I sources from this 
study. The error bar in the upper right corner shows the typical uncertainty in luminosity 
and effective temperature for a class I source. Uncertainties for T Tauri stars are ~ 33%-50% 
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of the uncertainty for class I sources. The arrow indicates the change in log L and log Te if a 

continuous source of emission veils the optical contimium. The dashed lines plot isochrones 
for the pre-main-sequence tracks of D'Antona & Mazzitelli (1994) at times of 10^, 10^, and 
10'' yr from top to bottom. The dot-dashed line is the stellar birthline from Stahler (1988). 
The solid lines are pre-main-sequence tracks for stars accreting from discs at 1O~^M0 yr~^ 
(thicker line) and 1O~^M0 yr~^ (thinner line). 
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Table 2. Emission line equivalent widths from FAST spectra 


Source Name 


StD 


r / r 


TT Ftvt ttI n r o a tt« t /^/^to To ttI n^-x^y i f"y^c\ 

Ha [N llj 6584 He 1 6678 [S llj 6717+6730 


HBC 351 


HI 


0.05 


2.8 


HBC 352 


HI 


<0.03 


-2.2 


HBC 353 


III 


<0.04 


—2.1 


HBC 354 


III 


0.14 


— 1.2 


TTTI O P* P* 

HBC 355 


HI 


0.07 


— 1.2 


HBC 356 


HI 


0.06 


0.8 


HBC 357 


HI 


0.10 


1.1 


HBC 358 


III 


0.04 


10.0 


TTT^ r» 

HBC 359 


III 


0.04 


5.4 


HBC 360 


III 


<0.20 


6.7 


TTT^ /~1 r» -1 

HBC 361 


HI 


<0.20 


6.4 


TTT~» o r» 

HBC 362 


HI 


<0.15 


6.6 


LkCa 1 


III 


<0.05 


3.4 


Anon 1 


III 


<0.02 


2.5 


04108+2910 


II 


<0.43 


37.5 ■ ■ ■ ■ ■ ■ 2.0 


V773 Tau 


II 


0.10 


2.8 


FM Tau 


II 


0.19 


101.0 ••• 1.6 


FN Tau 


II 


0.31 


16.0 


CW Tau 


II 


0.23 


-\ A r\ r\ -\ c\ -1-1 

140.0 ••• 1.2 1.1 


CIDA 1 


II 


<0.18 


149.0 


LkCa 3 


HI 


0.01 


2.7 


FP Tau 


III 


0.13 


27.0 0.7 


FO Tau 


II 


0.15 


116.0 ••• 1.5 


CX Tau 


II 


0.17 


18.0 


CIDA 2 


HI 


<0.09 


6.8 


LkCa 4 


HI 


<0.02 


3.2 


CY Tau 


II 


0.10 


63.0 ••• 1.7 


LkCa 5 


III 


0.16 


3.8 


CIDA 3 


II 


<0.14 


10.5 


V410 lau 


TTT 

111 


0.01 


2.3 


DD Tau 


II 


0.40 


177.5 ■ ■ ■ 3.9 2.2 


CZ Tau 


II 


0.29 


6.7 


HBC 372 


III 


<0.15 


— 1.2 


V892 Tau 


II 


0.58 


6.6 


Hubble 4 


HI 




3.9 


HBC 376 


HI 


<0.05 


1.9 


FQ Tau 


II 


0.15 


81.0 ••• 3.0 


BP Tau 


II 


0.10 


92.0 ••• 1.0 


V819 Tau 


HI 


0.03 


3.2 


LkCa 7 


HI 


<0.02 


3.9 


J2-157 


HI 


<0.50 


15.5 


DE Tau 


II 


0.11 


59.0 



Table 2. Continued 



Source Name 


oiljiJ 


r „ , „ / r , 
J^firI J^h 


TT/-./ 
IIQ 


[i\ llj 00o4 


rie i Do / o 


TQ TTl (^71 7-Lfi7^n 
[D llj Dili "tD I oU 


xti iau 


TT 


U.zo 


1 Q n 
lo.U 








rlJJ zooO ( z 


TTT 
ill 


no 
<U.Uz 


— l.o 








r o iau 


TT 

11 


U.oU 


iti n 
m.U 


O.O 


l.o 


Q n 

o.y 


i iau 


TT 

11 


n on 
U.zU 


/I n c 
4U.0 




1. ^ 


1 7 
1. / 


I'D Trin 

ir^ iau 


TT 
11 


U. ID 


1 n 
lU.o 








T1 /1Q70 
J l-4o ( Z 


TTT 
111 


<U.U4 


Q 

z.y 








r V iau 


TT 
11 


U.lb 


on n 
zU.U 








ur iau 


TT 
11 


n HQ 


04.0 




1. ( 




ukj iau 


TT 
11 


U.tJU 


1 o 1^ n 
Izo.U 


n 7 
U. / 


l.U 


Z.4 


iliJO ooo 


TTT 
ill 


n no 
U.Uz 


— l.D 








iJil iau 


TT 
11 


n 1 
U. Iz 


oo.O 




1 A 

1 .4 


U.D 


T\T Ton 

iJi iau 


TT 
11 


U.UtJ 


o n 
z.U 








1C\ Ton 

ii^ iau 


TT 
11 


n 1 o 
U. Iz 


7 7 
/ . < 








Ton 

rA iau 


TT 
11 


U.lb 


1 /I C 
14.0 








T\IC Ton 

iJiY iau 


TT 
11 


U. lo 


/1Q 

4y.O 




1 .u 


1 n 

l.U 


Ton 

lZj iau 


TTT 
ill 


U.U ( 


14. U 








TTJ 

Jil 00 


TTT 
111 


<U.UD 


o o 
z.z 








"\7"007 Ton 

vyz/ iau 


TTT 
111 


<U.Uo 


inn 
lU.U 








ili)V_j oaz 


TTT 
111 


U.Uo 


1 1 
1 . 1 








UT Ton 

ilij iau 


TT 

11 


U. ( ( 


/I Q n 
48. U 






o.D 


V7 Ton 

AZj iau 


TT 
11 


U.Ol 


1 7t^ n 
1 ( O.U 




1 .4 


o 
z.z 


UTZ Ton 

ilis. iau 


T /TT 
1/11 


n o7 
U.z 1 


OO.O 




U.4 




"\7''71 n Ton A 

V / iu iau A 


TT 
11 


U.Uo 


OO.O 








"\7''71 n Ton FJ 

V / iu iau i) 


TT 
11 


n 1 /I 
U. 14 


4.0 








ijiOOi-Oi 


TTT 
ill 


<U.U4 


1.0 








"V/OOQ Ton 

vyzo iau 


TTT 
ill 


<U.U4 


1 .o 








"\7"QO'7 Ton 

\oZi iau 


TTT 
111 


<U.U4 


A A 
4.4 








"\rQO^^ Ton 

V ozo iau 


TTT 
111 


n no 
U.Uz 


Q 1^ 
O.O 








ir"V Ton 

r 1 iau 


TT 

11 


n 1 o 
U. iz 


7Q n 
( o.U 




U.z 




XrV Ton 

r L iau 


TT 
11 


oc; 


1 Qn c: 
loU.O 




z.U 




Ton 

(jtj iau 


TT 
11 


n 1 o 
U. la 


ou.o 




n 7 




TTV Ton 17 

uZj iau iL 


TT 
11 


no 
<U.Uz 


DO.U 




U.D 




TTV Ton \HT 

uZj iau w 


TTT 
111 


no 
<U.Uz 


/I Q C 
4o.O 




U.D 




Jil iiz 


TT 
11 


n 01 
U.zl 


1 n 
10. u 








L1551-55 


III 


<0.06 


1.8 








GH Tau 


II 


0.08 


27.5 








V807 Tau 


II 


0.01 


13.5 








V830 Tau 


III 


<0.02 


1.8 








GI Tau 


II 


0.19 


14.5 




0.8 




GK Tau 


II 


0.13 


30.5 




0.9 




DL Tau 


II 


0.33 


101.0 




2.2 


0.4 


HN Tau 


II 


0.43 


145.0 


1.0 


1.7 


5.1 



Table 2. Continued 



Source Name 


Q1hT» 




IIQ 


rivr TTl f\^SiA 
[i\ llj O0o4 


TTq T (^(X7Si TQ TTl /^71 7j_^^7Qn 
He 1 DO to [O llj / 1 / "tO I OU 


1 1 lau 


TT 

ii 


<U.Uo 


O 7 






L^i lau 


TT 
ii 


U.l ^ 


7^^ 




1 n 
i.U 


JJM lau 


TT 
ii 


U.lD 


1 1 /I n 
1 14. U 




O 1 

z.l 


ITJ 1 no 

jn lUo 


TTT 


<<f] 1 1 


o.o 






±1 J_) Vy 4:U / 


TTT 
111 




1 .o 






AA lau 


ii 


U. lo 


Afi n 
40. U 




1 1 
1.1 


nu lau 


TT 
ii 


U.oo 


1 ni K 
lUl.o 




o o 
z.z 


r r lau 


TTT 
iii 


U.Uz 


z.o 






JJ1\ lau 


TT 
ii 


n no 

u.uy 


ion 
lo.U 






111^ lau 


TT 


U.OD 


on n 
zU.U 






m'o T-^ii /r^o 
Hr^ lau/Lrz 


TTT 
iii 


<U.Uo 


O Q 






Hr^ iau/L:ro 


TTT 
iii 


<U.U1 


O Q 
Z.O 






waro o-zo 


T 
i 


U.oz 


01.0 






T VPn 1 zL 


TTT 

111 


U.Uo 


1 1 
1 . i 






r»n T;:m 


TT 
11 


n 9zL 
u.z^ 


1 01 n 

±U± .U 




9 14 

Z.U 


ITV Tan 


TTT 
111 


\ U.UZ 


o . u 






V 1 lau 


TTT 
111 


n 1 /I 
U. 14 


7 Q 
/ .O 






_LKL>a lo 


TT 
11 


nil 
U.ll 


lo.O 






TTT OOQ 

Jll ZZO 


TT 
11 


<U.l i 


/I O 

4.Z 






TW Ton 

1 vv 1 au 


TTT 
111 


no 
•CU.Uz 


A n 
4.U 






L-oivu iau/4 


TT 
11 


U.oo 


O Q 

z.o 






ijKnaooz / KjZ 


TTT 
ill 


no 
<U.Uz 


O.O 






IjKllaoOZ/Ljl 


TT 
11 


U.Uo 


4.y 






vyoo lau 


TT 
11 


U.lO 


on n 
zU.U 






Ur lau 


TT 
11 


U.ZD 


1 ni K 
lUl.o 


1 1 
1.1 


1.4 4.0 


LtU lau 


TT 
11 


U. 1 < 


TQ n 
t o.U 




1 Q 
i.O 


\Jy^ lau 


TT 
11 


U.zl 


1 no n 
lUz.U 






riaro u-o( 


TT 
11 


nil 
U.ll 


1 Q n 
lo.U 






T\T> Ton 

UK lau 


TT 
11 


n OQ 
U.zo 


7Q n 
/o.U 




1 K 
1.0 


lau 


TT 
11 


n 1 n 
U. lU 


OO.O 




1 1 
1.1 


TTV An*. 

u I Aur 


TT 
11 


U.oo 


Oo.o 




no 1 Q 
u.y i.O 


ot o4 


TT 
11 


U.lO 


on A 

yu.4 




O A 

Z.4 


LriVi Aur 


TT 
11 


n on 
U.zU 


/ y.o 








TTT 
111 


<rf\ 01 

\ U .U 1 


1 9 

1 . z 






AB Aur 


II 


<0.19 


44.0 




0.5 


SU Aur 


II 


0.18 


2.0 






HBC 427 


III 


0.01 


1.4 






V836 Tau 


III 


0.12 


7.7 






CIDA 8 


III 


<0.11 


79.0 




2.7 


CIDA 11 


III 


<0.13 


49.0 






RW Aur 


II 


0.18 


75.5 




1.2 1.4 



